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Abstract 

The spin-diffusion constant of the 2D t — J model is calculated for the 
first time using an analytical approach at high temperatures and a recently- 
developed numerical method based on the Lanczos technique combined with 
random sampling in the intermediate temperature regime. A simple relation, 
a = DsXi between spin conductivity and spin diffusion is established and 
used to calculate the latter. In the high-temperature and low-doping limit 
the calculated diffusion constant agrees with known results for the Heisen- 
berg model. At small hole doping, Dg increases approximately linearly with 
doping, which leads us to an important conclusion that hopping processes en- 
hance spin diffusion at high temperatures. At modest hole doping, 5 ~ 0.25, 
diffusion exhibits a nonmonotonic temperature dependence, which indicates 
anomalous spin dynamics at small frequencies. 
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I. INTRODUCTION 



There appears to be a rising interest in the spin dynamics of high-temperature supercon- 
ductor compounds. The research in this field is mostly concentrated around the spin-lattice- 
relaxation rate, 1/Ti. In particular, recent measurements of 1/Ti in Sr-doped La2Cu04 
display rather unusual temperature and doping dependences. Sokol et al. have calcu- 
lated 1/Ti for the undoped case and found good agreement with experiments. They also 
calculated the spin-diffusion constant, which gives a short-wavelength contribution to the 
nuclear-relaxation rate at high temperatures. The relative strength of this diffusive contri- 
bution to the relaxation rate depends strongly on the hyperfine structure, which differs in 
the case of Cu and O nuclei Spin diffusion can serve also as a general numerical test 

probe for the spin dynamics of the system at small frequencies and momenta. 

The spin-diffusion constant of a paramagnet in the limit of high temperature was first 
calculated by de Gennes using a moment expansion of the spin-susceptibility function 0. 
The method was later improved by Bennett and Martin using an exact spectral repre- 
sentation of the spin-correlation function in the hydrodynamic limit. Morita investigated 
a theoretical criterion for the occurrence of spin-diffusion by introducing a friction function. 
He also calculated diffusion constants for various lattices and values of spin. 

In this paper we present a finite-temperature calculation of the spin-diffusion constant 
in the t — J model on a square lattice 

H = -t J: (4cj, + h.c.) + J ^ (Si ■ Sj - n,n^/A) , (1.1) 

<ij>,s <ij> 

where < i,j > represents summation over pairs of nearest neighbors, c\g and are cre- 
ation and annihilation fermion operators restricted to the basis of single-fermion occupation 
number, and Si is a spin-1/2 operator. 

To obtain finite-temperature static and dynamic properties of the t — J model, we use a 
recently developed numerical method, based on the Lanczos exact-diagonalization technique 
combined with random sampling First we present an alternative method of calculating 
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spin diffusion, which is based on the spin-conductivity function. Next, we present an analyt- 
ical high-temperature result for Dg in two limits of the t — J model: a) t = and b) J = 0. 
In the case of zero hole-doping, 5 = 0, our analytical method reproduces the known result 
for Ds of the Heisenberg model. Next, we briefly discuss the finite-temperature Lanczos 
method, and then present results for Dg in the t — J model at arbitrary parameter values. 
Last, we present conclusions. 



II. THE METHOD 

A. Spin-difFusion constant 

We start by noting that the total 2;-component of the spin operator = J2rS^{r) is a 
constant of motion from which follows the operator-continuity relation 

^^^ + V-j.(r,t) = 0, (2.1) 

where is the spin-z-current operator. Now, we suppose that at finite temperatures the 
hydrodynamic description can be applied to spin relaxation in the t — J model in which the 
net flow of magnetization from the region of large towards small is governed by the 
diffusion equation 



d{S {r,t))non-eq _ p. ^2 I oz 

-i-^s V \D yi- 1 J / non—eq 



DsV^{S^{r,t))r,on-eq = 0. (2.2) 



It is important to notice that contrary to eq. ( |2.lD , which is an operator equation, the 
diffusion equation is valid only with brackets Onon-eq, which represent nonequilibrium 
ensemble averages. 

Applying linear response, hydrodynamics and the fluctuation-dissipation theorem one 
can obtain the following form of the Fourier transform of the spin-correlation function 
S{r,t) = {S^{r,t)S%0,0)} [|], 

S{k,uj) = {S'{k,u)S'{-k,-u)) = 

2 ujDsk'^x / x 

~ X (2 ^) 

l_e-/3- u^ + {Dsk^y' ^ ■ ' 
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valid for small k and u. The 5'^(k, u) is a spacial and temporal Fourier transform of S^{r, t). 
We use units in which h = = 1- In equation ( |2.3| ) (3 is the inverse temperature, x is 
the static magnetic susceptibility, and Dg is the spin diffusion constant. Note, that the 
last factor in eq. ( |2.3| ) is the imaginary part of the dynamic susceptibility function in the 
hydrodynamic limit. The brackets <> in equation ( |2.3| ) and all the following equations 
represent thermodynamic equilibrium averages. According to Onsager, the spontaneous 
equilibrium fluctuations of the spin, described by the spin-correlation function ^(k, u), obey 
the same diffusion equation as do the nonequilibrium - induced spin fluctuations in eq. ( p. 21) . 
Next, we define the spin conductivity as 

a^(k,uj) = — {j^{k,uj)j^{-k,uj)) = 

2u 



1 - e-^'^ 



/oo 
dtJ2e'i^---'\Ur,t)U0,0)), (2.4) 
-oo J. 



2uj 

where jx{r,t) is the x— component of the spin current and j^(k, cj) its Fourier transform. 
The spin conductivity o"(k, uj) can be related to the spin response function ^(k, uj) using the 
continuity relation in (k, cu) space: —ujS''(k,Lj) + k ■ j(k, cj) = 0: 

a^{k,u) = luil-e-^n^^- (2.5) 

Combining eqs. ( p. 3] ) and ( |2.5| ) we derive the following simple Einstein relation between the 
spin conductivity and spin-diffusion constant: 

a = a4k = 0,cu = 0) =D,X- (2-6) 



B. Dg in the high temperature Umit 

In the high-temperature limit of the Heisenberg model it is reasonable to expect that the 
spin conductivity, cr(cij) at k = 0, can be approximated with a Gaussian form 0. We have 
checked this assumption numerically. However, this is in general no longer true in the case of 
the t — J model, predominantly due to the coexistence of two different energy scales. There 
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are two limiting cases of the t — J model, however, for which we expect that a Gaussian 
form is still a reasonably good ansatz: a) J = 1, t = 0, 5 <^ 1 and b) J = 0, t = 1, 5^1, 
where 5 represents hole-doping. Note, that the above assumption must be limited to small 
hole-doping regime. In the high-doping regime, the spin conductivity has a Drude-like 
(Lorentzian) form. 

While a direct computation oi a{uj = 0) would be subject to all the usual vagaries of 
numerical computation, we can avoid these (after taking the assumption of the Gaussian 
form of <j{uj)) by taking a ratio (see ( p.lO|) ) of the first two nonzero moments of the spin 
conductivity: 



/OO jr 
du;a{uj) = -J2{Ur,0)U0,0)) 
OO 1 J. 



^d.;.;M^) = -E((^^)j.(r,t)(-^^)j.(0,t))), (2.7) 

where the time derivatives are to be taken at t = and the index r runs over all lattice sites. 
The spin-current operator ^^^(r) for the t — J model can be calculated using the continuity 
relation ( |2.1| ) together with the relation ■^S^{r,t) = —i[S^{r,t),H]. Taking into account 



the fact that all the operators are now defined on a discrete lattice, we obtain 

J [S^S;^^ - h.C.) - t (4^Cr+^-f - 4^Cr+5i| - h.c.)] , (2.8) 



t r 



with X being the unit vector in the x-direction. To avoid confusion we note again that 
represents the x-component of the spin-^-current operator. In order to calculate the first two 
moments ( |2.7] ), one has to compute also the time derivative of the spin current. Taking the 
high-temperature expectation values we arrive after a somewhat tedious but straightforward 
calculation at the following expressions for the first two moments: 



(-°)=(^^(l-5) + y(l-5)j ^, 



{uj') =i^^6{l-6f{9 + 6) + ^{l- 5f (5 - 35) + 0(tV^)) ^. (2.9) 
In the above equation, the second moment {uo"^) is correct only in the two limits of interest. 



i.e. at a) J = or b) t = 0. Assuming a Gaussian functional form for the spin conductivity, 
we can express the diffusion constant in eq. ( p^ ) in terms of the first two moments ( p.7| ) 



1 1 



^ ^ (2.10) 



which leads to the final expressions for the spin-diffusion constant for the two limiting cases: 

D,= { ' (2.11) 

(9+^x1-^) ; = 0- 

To test the method, we compare the diffusion constant (|2.11| ) in the Heisenberg limit (t = 
0, 5 = 0) with the known result for the Heisenberg model at high-temperatures Dg/ J = 

and find exact agreement. At t = and small hole-doping 6, Dg decreases linearly with 
doping as Dg/J ~ J^i^ ~ "m'^)- This is simply related to the fact that the number of 



20V 10' 

spin-current carriers (spins) decreases with doping. In the opposite limit, when J = 0, spin- 
diffusion ( p. lip increases linearly with doping as Dg/t ~ (2vT0/3)y^(5. It is straightforward 
to check that at zero doping and J = 0, Dg/t = 0, since the spin-current ( p.8| ) at 5 = is 
zero. The linear increase of Dg with doping S can in this case be explained by the linear 
increase of spin-current carriers (holes). Interestingly, the coefficient of the linear increase 
of Dg in the J = limit is larger than the coefficient of decrease in the t = limit. Since 
the physical parameter region (relevant for cuprates), where J/t ~ 0.3, is not accessible by 
the present analytic method, we make the following assumption: neglecting the combined 
hopping and spin-flip contributions to the diffusion constant we predict that in the high- 
temperature regime and at J/t = 0.3, Dg increases linearly with doping. 



C. Numerical method 

Having exhausted the possibilities of analytical calculation we now turn to the numerical 
calculation of Dg in physically relevant parameter regimes. We used a finite-temperature 
Lanczos method, recently developed by one of the present authors (J.J.) in cooperation 
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with P. Prelovsek. Since the details of the method have been pubhshed we only briefly 
summarize the main steps of the method. 

Taking into account the translational symmetry of the system, the spin conductivity 
at k = can be written as 

= / e-^-*(J,.(t)J,(0)), (2.12) 



2Nuj 

r 

where is the number of sites. According to the finite-temperature Lanczos method [Q, 
the thermodynamic equilibrium average in eq. ( p. 121) is expressed in terms of functions 



l^m)' iV'fc) and corresponding energies e„m, Cnfc, generated by the Lanczos method on a finite 
cluster with sites, 

No M 

{Ut)JM) = ^-^EE(^IOe-^^"'"e'^"'"''"^^*(d^xlC)(Cl^xl^), (2.13) 

" rn,k 
No M 

n m 

where states \n) for n = 1 to Ngt represent a complete orthonormal basis set of the Hamil- 
tonian H. In the first summation, only A^q ^ ^st of states \n) are randomly chosen for 
partial random sampling of basis states. The functions li'l^) and energies e^m are obtained 
as eigenvectors and eigenvalues of the tridiagonal form of the Hamiltonian, generated by the 
Lanczos procedure, 

H\rj = bnm\rm-l) + ^^nmlrj + 1 01^+1 ) (2.14) 

m = 0, . . . , M < Nst; &„o = &n(A/+i) = 0, 
starting from the initial wavefunction |0q) = \n). Similarly, {ip^) and ink are generated 



by the Lanczos procedure, analogous to (|2.14|), but with a different initial wavefunction 



m = J.\n)/^{n\J^Jn). 

The static magnetic susceptibility x = l3S(k = 0) is evaluated with the same method, 
adapted to the calculation of static quantities. Here, the equal-time spin correlation function 
is expressed as 



iVo M 

5(k) = ^-^EEe-^"'"(^IO(CI^^(k)5^(-k)|n), (2.15) 

n m 

where S^(k) = N~^^'^ X^r ^ '"'^'^(i"), and the rest of notation is the same as in Eq. (|2.13| ). 

The method yields very accurate results even for a severely reduced number of Lanczos 
steps M <^ Nst and number of random states A^^o <^ Nst- Provided that No is large enough, 
ax{uj) has correct frequency moments (u^) in the limit T ^ oo for p = 0, M (for a finite 
cluster). Similarly, at finite temperatures the corresponding double series for the moments 
of cr(c(j) are correct up to /?*'+^(c<j^^^^), for p = 0, . . . , M. The method also gives reliable 
results at T = 0, assuming that the number of Lanczos steps M is larger than the number of 
Lanczos steps required to obtain the corresponding ground-state results. For a typical run, 
M = 100 — 150 Lanzos steps and sampling over A^^o = 200 — 300 were sufficient to obtain 
reasonably accurate (within 5%) results for spin conductivity and static spin susceptibility. 



III. NUMERICAL RESULTS AND CONCLUSIONS 

Using the finite-temperature Lanczos method we computed the spin-diffusion constant 
of the t — J model in the physically relevant parameter regime, i.e. at J/t = 0.3 and for 
four different hole-dopings, 6 = 0, 1/16, 1/8 and 1/4. The calculation was performed on a 
4x4 cluster with periodic boundary conditions. 

In Fig. 1 we present numerical results for the renormalized spin-diffusion constant Dg/J 
of the Heisenberg model as a function ofT/J. For comparison we present on the same graph 
Dg/J obtained by the high-temperature expansion by Sokol et al, Dg/J = y^(l + 
with j3 = J/T. Numerical results agree well with the high-temperature predictions down 
to temperatures T/J ~ 1.5, where Dg obtained by the numerical method saturates. In 
the T — > limit, Dg has to diverge, since the Heisenberg model at T = has finite. 



nonzero spin stiffness [^. Saturation of Dg is therefore a finite-size effect. Nevertheless, the 
agreement with the high-temperature expansion at higher temperatures is quite encouraging. 
In contrast to the finite-size T = calculations, where results are strongly size-dependent. 
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the finite-temperature Lanczos method even on a small cluster size provides results which 
are correct in the thermodynamic limit. We expect results to be size independent as long as 
the temperature is larger than the average discrete level spacing e, which is a consequence 
of the finite-size system. At finite doping, where the low-lying energy sector becomes more 
dense, we expect size-independent results down to much smaller temperatures. 

We now turn to a doped case. In Fig. 2 we plot Dg/ J as a function of temperature for four 
different values of doping at J/t = 0.3. As predicted by the infinite-temperature moment 
expansion, Dg/ J is at high temperatures, T ^ t, increasing approximately linearly with 
doping 6 foT 6 <^ 1. In contrast to the undoped case where Dg as a consequence of finite-size 
effects saturates around T J, Dg at finite doping diverges with decreasing temperature 
(even in the finite cluster). This is an indication of the formation of a quantum coherent 
state, which has a finite spin-stiffness at T = 0. Dg at small doping, 6 < 1/8, displays 
a monotonic l/T-like temperature dependence, similar to the undoped case. Somewhat 
unexpected is the nonmonotonic temperature dependence of Dg at moderate doping, 6 = 
1/4. At high temperatures, Dg/J decreases with falling temperature, reaches a minimum 
around T ~ 1.5 J, and when the quantum coherent regime sets in at even lower temperatures, 
sharply increases. We attribute this unusual nonmonotonic temperature dependence of the 
diffusion constant to the increase with doping of the number of low-energy short-lived many- 
body states. An enhanced density of low-energy states contributes to enhanced inelastic 
scattering, which consequently leads to a decrease of the diffusion constant. This picture is 
in agreement with the enhanced entropy at low temperatures |Tl| and strong charge-carrier 
scattering, which is manifested in the linear temperature dependence of the resistivity in the 



t - J model 12 



In conclusion, we have calculated the spin-diffusion constant Dg for the t — J model. In 
the high-temperature limit, Dg can be approximately calculated via the moment expansion 
of the spin-conductivity function at t = or J = and for small doping 6. In the Heisenberg 
limit, this method gives the correct result for Dg. At small doping, 5 <^ 1, and at t = (the 
dilute Heisenberg model), Dg/J decreases linearly with doping, indicating that the spin- 



current carriers are localized spins. At J = and 5^1 Ds/t is increasing linearly with 5 
which is consistent with the fact that the spin-current carriers in this limit are mobile holes. 

At finite temperatures, we used a recently developed Lanczos method to obtain Ds as 
a function of temperature and doping. We found good agreement of Ds in the Heisenberg 
limit with the high-temperature expansion result at higher temperatures. In agreement with 
our analytical predictions, Ds is increasing approximately linearly with doping at Jjt ^ 0.3 
in the high-temperature and low-doping regime. Since the diffusion constant decreases with 
hole doping in a dilute Heisenberg model, i.e. at t = 0, we can conclude that hopping pro- 
cesses contribute substantially to spin diffusion at high temperatures. At moderate doping, 
5 = 1/4, and high temperature, T ~ J, decreases with decreasing temperature. This 
anomalous behavior can be associated with the recently reported anomalous spin dynamics 
of the doped t — J model at low-frequencies and enhanced entropy at low temperatures 
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FIGURES 

FIG. 1. Spin-diffusion constant Ds of the Heisenberg model vs. temperature obtained 
by the finite-temperature Lanczos method (full circles). For comparison we present the 
high-temperature-expansion results obtained by Sokol et al. (dashed line). 

FIG. 2. Spin-diffusion constant Dg of the t — J model vs. T/t at J/t = 0.3 and various doping 
values. The arrow above the temperature axis points to T = J. The accuracy of results at finite 
doping is 5%. 
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